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I. INTRODUCTION 



Solving the Schrodinger equation with time-dependent 
boundary conditions, including moving ones, is a very 
hard work and can only be done in a few cases. See, e.g., 
[iHl] and Q for a recent review. Very interesting effects 
arc seen in such problems; diffraction in time [5[ is just 
such an instance. This phenomenon is characterized by 
quantum temporal oscillations in matter waves released 
from a confining region. It was shown when a wall acting 
as a perfect mirror, moves with finite velocity along the 
direction of propagation of a beam, the visibility of the 
fringes is enhanced . Moshinsky's theoretical work has 
been extended to the case of particles which are suddenly 
released from a ID box 0] and to particles with angu- 
lar momentum |8j] initially trapped in a hard spherical 
box. Exact solutions of the Schrodinger equation for a 
particle in a ID box with a moving wall have been found 
[3, • Using the semiclassical approximation, Luz and 
Cheng evaluated the exact propagator of the prob- 
lem. Grosche [llj did this task independently by means 
of an exact (mid-point) summation of the perturbation 
series with point-like perturbations. 

The motivation of the present work is that sudden re- 
moval of the boundary is an idealized case, thus we con- 
sider a moving boundary instead of a sudden removal of 
the boundary. The limit of infinite velocity of the mov- 
ing boundary clearly corresponds to the sudden removal 
case. Possible applications to optical effects connected 
with moving mirrors are additionally encouraging. For 
these reasons we aim to solve the Schrodinger equation 
for a particle in a hard sphere with varying radius. 



II. EXACT SOLUTION 



Consider a particle with mass fi inside a hard sphere 
with a time-dependent radius L(t). The potential energy 
function is zero if r < L{t) and infinite otherwise. The 
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Schrodinger equation is then 



(1) 



with the boundary condition 'ii{r,t)\r=L{t) = 0. 

The instantaneous energy eigenfunctions and eigenval- 
ues are respectively 
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I = 0,1,2,...; n ^ 1,2,3,... and m = -1,-1 + - 
1,1, where ji{x) and Yim{9,(j)) are respectively spherical 
Bessel functions and harmonics, xin is the ri*'^ zero of the 
spherical Bessel function of order /, i.e., ji{xin) — Q. It 
must be noted that all Bessel functions with I ^ have a 
zero at the origin, but to have a non-zero wave function 
these zeros must be excluded. 

Using the method of "separation of variables" for solv- 
ing the partial differential equation ([T]), we propose the 
solution 
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where we have used the spherical symmetry of the Hamil- 
tonian. 

Putting eq. (U) into eq. ([1]) one gets 
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The radial part of the proposed wave-function, R{r, t) — 
U{r,t)/r, must be zero on the shell, thus the boundary 
conditions on U{r,t) are U{r,t)\r=o = = U{r,t)\r=L{t)- 
Now, we follow 01 to solve the eq. By defining a 
new coordinate 



we get 
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where L(t) = dL{t)/dt and moving boundary conditions 
are replaced by fixed-boundary ones; U{s,t)\s^Q = = 
U{s,t)\s=i- When the transformation 
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is introduced in eq. (O, one obtains 



2fiL^(t) [ ds 



^{s,t) [9) 



for the uniform motion of the wall, i.e., L{t) = 0. Bound- 
ary conditions on if{s, t) are if{s, t)|s=o = = if{s, t)\s=i. 
Defining the new time variable r as 
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Eq. © transforms to 
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Inserting (p{s,T) = exp{~iE'T/h)il;{s) in ([TT|) . one gets 
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By introducing new variable k'^ = 2iJ,E' /H^ , we obtain 
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The solutions of this equation arc spherical Bcssel func- 
tions 



■)p{s) = s[ciji{ks) + C2ni{ks)] 
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If the radial wave-function R{r) is finite at the origin, 
C2 = 0. The requirement that ip^s) = at s = 1 means 
that k can take on only those special values 

kin = xin (n = 1,2,3, ...) ■ (15) 

For the uniform change of the radius with velocity u 

L{t) = a + ut , (16) 

where a is the initial radius, one has 
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By using equations (H?]), ((11]), ^ and © one 

obtains 
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for the radial part of the wave-function. Unknown coef- 
ficient Ci is determined by the normalization condition 



J drr"" J dn\-^inm{r,t)\^ = 1 , 



where 



'finmir,t) = Rin{r,t)Yi„,iO,(f>) 



(19) 
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are the solutions of the Schrodinger equation ([T]) for a 
particle in a spherical box with a wall in uniform motion 
and J dfl ^ J_-^d{cos9) J^^^ d(j). 

Using the orthogonality of the spherical Bessel func- 
tions |9l 
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one obtains 
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Thus apart from a phase factor, one obtains 
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where we have introduced new dimensionless parameters 
a = iiau/{2h) and ^(t) = L{t)/a. 

Functions \I'i„m(r,i) vanish at r = L(t), remain nor- 
malized as the radius changes, and form a complete or- 
thogonal set. The general solution of eq. ([T|) is a super- 
position of functions 



oo oo 



*(r,t) = XI H H Cl>n'm''i'l'n',n'{r,t) , (24) 



I'—O 7i' — l ■m' — — l' 



with time-independent coefficients ci'n'm' determined 
from the relation 
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General solution can also be expanded in terms of in- 
stantaneous eigcnfunctions as 
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now with time-dcpcndcnt cocfRcients biin'm'(t) deter- 
mined from the relation 
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Using eqs. p7p and ([24| and the orthogonality of spher- 
ical harmonics, one finds 
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This integral is not elementary and following the proce- 
dure of [Sj, can be reduced to a combination of terms 
involving the Ferensel integrals and derivative of Legcn- 
dre polynomials. 

The expectation value of the energy of the particle is 
obtained from 

{E{t)) = \bl'u'm'it)\^ Ei,n-m'[t) . (30) 

V n'm' 

If the particle is initially in an energy eigenstate, i.e., 
^'(r,0) = u,„,„(r,0), then 
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which is not an unexpected result as quantum numbers 
I and m do not change. 



III. NUMERICAL CALCULATIONS 

Numerical computations are shown in figs. [T] and [5] for 
a particle that is initially in the first excited state with 
three fold degeneracy. In this case we have 
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FIG. 1: (Color online) Transition probabilities versus £^{t) for 
different values of velocity parameter a: a) a = —2; h) a = 
—4; c) a = —6 and d) a = —10. In each part the black curve 
shows |6iimP, red one |fei2mp, green one |6i3mP and the blue 
one 6i4mP. 




FIG. 2: (Color online) Ratio of the energy expectation value 
to the instantaneous first excited energy as a function of ^(t) 
for three different values of velocity parameter. 



for the ratio of energy expectation value to the instanta- 
neous first excited state energy. 

Figure [T] shows the squares of energy eigenfunction ex- 
pansion coefficients versus ^(t) for three different con- 
traction rates a. For these values of a, it was found that 
series (pS)) converges for the first ten terms. 

Figure [5] shows the ratio of the expectation value of 
the energy to the energy the particle would have if it 
remained in the first excited state Unm for the sphere 
in contraction. Here fifteen terms in eq. p2p leads to 
convergency. 



We have plotted dimensionless radial probability den- 
sity piniviniTin) = Af„r;f„ |i?(77;„ , T/„) | ^ in fig. El for a 
particle initially in the state uo,5.o, against dimensionless 
position coordinate rjin ~ r/Xin at dimensionless time 
coordinate Tj^^-* = i^;„t*-°-' and in fig. 0] for a particle ini- 
tially in the state (a) Mo.15,o and (b) uo,ioo.Oi against di- 
mensionless time coordinate Tin = i^int at dimensionless 
observation point Ty^^' = r'^^^/Ajn; where Az„ = 27ra/a;;„ 

and i^in — Ein/h. In our calculations r'*'^ = 2a and 
i(0) ^ (7.(0) 
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FIG. 3; (Color online) Dimensionless radial probability den- 
sity pin(jlin,Ti„) for a particle initially in the state ?io,5,o, 
against dimensionless position coordinate ■qin at dimension- 
less time coordinate t/^' , for six different values of expansion 
rate; a = Q (black curve), a = O.Olain (red curve), a = ain 
(green curve), a = 10ai„ (blue curve), a = 15q:i„ (yellow 
curve) and a — 20ai„ (magenta curve); where ain = xin/2. 



FIG. 4: (Color online) Dimensionless radial probability den- 
sity pinivintTin) for a particle initially in the state (a) 
^^0,15,0 and (b) mo, 100,0, against dimensionless time coordi- 
nate Tin = J^int at dimensionless observation point ?7;^', for 
three different values of velocity parameter; a = 0.9ai„ (black 
curve), Q — ain (red curve), a — 2ai„ (green curve); where 
ctin ~ xin/2. Ti and T2 are dimensionless classical flight times 
from the front and back edges of the sphere to the dimension- 
less observation point. 



Ti and T2 arc dimensionless classical flight times from 
the front and back edges of the sphere to the dimension- 
less observation point rj^^^ for a particle in the state uinm- 



where we have introduced the propagator as, 



IV. PROPAGATOR 

One can construct the propagator as follows 
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where S{t, Iq) is the time evolution operator and we have 
used the fact that if the particle is in the sate |^'/„m) 
at to, it remains in that state as the wall moves, i.e., 
S(t,to)\'^inmito)) = |*;rmi(t))- Now, wc writc this equa- 
tion in the form 



One sees when I = 0, eq. ([7]) reduces to eq. (4) of (H, 
i.e., I = corresponds to a particle in a ID box with the 
left wall at X = and the right wall in uniform motion. 
In order to have the relation 
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in ID, we must write ID propagator as 



SUMMARY AND DISCUSSION 
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where we have used Yqq ~ l/v47r. The first two Bcssel 
functions are 
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tlius .Ton = mr and jiixon) — (— 1)"+ /jitt. Using these 
in eq. ([37]), we find 
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Now from eq. ([55)) we obtain 
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which is exactly eq. (32) of rcf. [l^ for the propagator 
of a particle in a ID box. This equation can be written 
in a compact form as a combination of 1^3 functions 



In this letter we found solutions of the Schrodinger 
equation for a particle confined in a hard spherical 
trap with a moving wall at constant velocity. We 
see in solutions ([25]) . except for the phase factor 
exp (— I / dtEinmit)/h) which has no coordinate depen- 



2ri L(t) 



ap- 



dence, a coordinate-dependent phase exp 

pears. It has been shown that this factor leads to an 
effective quantum non-local interaction with the bound- 
ary: even though the particle is nowhere near the walls, 
it will be affected 

Ftom fig. [TJ one sees that as the velocity of the wall 
increases, larger amounts of energy states other than the 
initial one, i.e., uum, are mixed in. Fig. [5] shows that 
for rapid contraction, energy expectation value increases 
faster than the 1/L^(t) increase which would be obtained 
in a quasistatic contraction. These results arc in agree- 
ment with the ones of ref. 0] obtained for a particle in 
an infinite square well with one wall in uniform motion. 
Confinement of the particle to a smaller region leads to 
enhancement of the energy expectation value. This can 
be explained by an application of the "old quantum the- 
ory" [IJl or by uncertainty relations |14|. 

In the process of expansion, there are two characteris- 
tic times involved: te, over which the parameters of the 
system change appreciably, and ti, representing the mo- 
tion of the system itself. In our calculations, tc = a/u and 
ti = a/vin. Figure [3] shows that for tc ^ (u ^ «;„), 
the particle, initially in the state ^0,5,0, will end up in 
the corresponding state of the expanded well. This pro- 
cess characterizes an adiabatic one for which external 
conditions change gradually [l^. While, in the oppo- 
site limit, rapidly changing conditions prevent the sys- 
tem from adapting its configuration during the process, 
hence the probability density remains almost unchanged. 

Noticing fig. SI one sees a quasi-classical behavior in 
the high-energy limit [l6j as the velocity of the wall in- 
creases. A non-monotonous increasing behavior of the 
density is seen for T < Ti only when u > vin, while for 
T > T2 a non-monotonous decreasing behavior is seen 
irrespective of the wall velocity. These results are in con- 
trast to classical mechanics. The height of the first max- 
imums decrease with u. The constructive interference 
with the reflected components from the wall for u < vin 
leads to this enhancement. Long time behavior of the 
density in the given observation point, is the same for all 
values of the wall velocity, which is not an unexpected 
result noticing the behavior of functions "^inm at long 
times. 

Propagator of the problem was derived using the 
spectral decomposition. 
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